Spin transport in diffusive superconductors 
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We employ the Keldysh formalism in the quasiclassical approximation to study transport in a 
diffusive superconductor. The resulting 4x4 transport equations describe the flow of charge and 
energy as well as the corresponding flow of spin and spin energy. Spin-flip scattering due to magnetic 
impurities is included. We find that the spin-flip length is renormalized in the superconducting case 
and propose an experimental system to measure the spin-accumulation in a superconductor. 
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Manipulation of spin-polarized currents can be used 
to study fundamental transport processes and might also 
provide new functionality in electronic devices. In fer- 
romagnets (F), the current is spin-polarized due to the 
spin-dependent density of states and the spin-dependent 
scattering potentials. In contrast, in s-wave supercon- 
ductors (S), electrons with spin up and spin down and 
opposite momentum form Cooper pairs with no net spin. 
Nano-scale superconductors therefore display strikingly 
different properties when driven out of equilibrium by 
spin transport than by charge transport. 

Most of the recent activities on the transport proper- 
ties of F/S junctions have studied effects caused by the 
physical properties on the F side of the junction. The 
zero spin Cooper pairs prevent spin-polarized electrons to 
flow into S. Consequently, a spin polarized current from F 
injected into S can result in noncquilibrium spin accumu- 
lation near the F/S interface. The competition between 
electron-hole correlations and spin accumulation on the F 
side has recently attracted considerable interest [1]. Pos- 
sible influence of the ferromagnetic order parameter on 
the superconductor has received less attention. Singlet 
pairing does not allow a spin accumulation in the super- 
conductor. Consequently, spin accumulation can reduce 
the superconducting gap and change the transport prop- 
erties both for transport via quasiparticles and for the 
supcrcurrcnt. Experimentally, spin transport in diffu- 
sive S has recently been studied [2]. Here, the reduced 
quasiparticle penetration due to spin accumulation re- 
sults in loss of spin memory which can be measured as a 
decreased magnctorcsistancc. 

Although the theory of nonequilibrium superconduc- 
tivity is widely used and developed, it has to the best of 
our knowledge not been completely generalized to study 
spin transport. In this work we thus use the Keldysh for- 
malism and the quasiclassical approximation [3-5] to rig- 
orously obtain a set of equations describing the transport 
of charge and energy in a diffusive weak coupling S, as 
well as the transport of spin. This will describe the pen- 
etration of spins into S and the associated suppression 
of the superconducting order parameter. Our descrip- 
tion of the transport properties will be based on a 4x4 
matrix equation formalism to include spin accumulation 



as well as electron-hole correlations. Spin-flip scattering 
from magnetic impurities is included as the dominant 
spin relaxation process inside the superconductor. We 
find that the spin-flip length is renormalized in the BCS 
case, and propose an experimental system to measure 
the properties resulting from the superconducting corre- 
lations. Many, but not all, experimental systems involv- 
ing spin-transport in superconductors are in the elastic 
transport regime [6], which is considered here. Comple- 
mentary studies based on the Boltzmann equation for 
spin-transport by quasiparticles in the inelastic trans- 
port regime have recently been published [7]. Note that 
spin-injection is qualitatively different in these opposite 
transport regimes due to the strong energy-dependence 
of quasiparticle flow in superconductors [7] . 

Let us now outline the derivation of our main results. 
We use natural units so that h = fee = 1, and the electron 
charge is e = — |e|. To describe the out-of-equilibrium 
electron-hole correlations as well as spin accumulation, 
we define the Keldysh Green's function as 

G«(l, 2) = ]T(-i) (fo) ik ([(V(l)) fc , (^(2)),-]-) , (1) 

k 

where ip = bPh V\b V 1 }' ^|] T * s a f° ur ~ vec t° r and i/j^ 
the corresponding adjoint vector. The matrix is 
the third Pauli matrix generalized to 4x4 space, p3 = 
diag(l, 1, — 1, — 1). The coordinates are 1 = (ri,ti) and 
2 = (r 2 ,i2). Similarly, we define 4x4 retarded and ad- 
vanced Green's functions (G R , G A ) in spin- and particle- 
hole space. 4x4 matrices are denoted by a "hat" super- 
script. A compact notation can be obtained by construc- 
tion of an 8x8 matrix in the Keldysh space (denoted by 
a "check" superscript) [5]. 

The quasiclassical Green's function is defined by 
g(R,T,p F ,E) = i /d£pG(R.,T,p,£T). This function 
is determined by the Eilenberger equation which in the 
mixed representation for a stationary state can be writ- 
ten 
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Here d = VI — ieA/53 is the gauge invariant derivative, 1 
is the 4x4 unit matrix, <f> is the electromagnetic scalar po- 
tential, A contains the superconducting gap and a is the 



2 



sclf-cncrgy due to clastic impurity scattering and spin- 
flip scattering by magnetic impurities in quasiclassical 
approximation. In the case of strong impurity scattering 
(dirty limit) transport is diffusive. Expansion of the qua- 
siclassical Green's function in spherical harmonics then 
gives the Usadel equations. The symmetries and normal- 
ization of the Green's function allows for a parameteri- 
zation of the quasiclassical, retarded component [3] 



1 cosh(6>) 
if2 sinh(6*)e _lx 



1T2 sinh(0)e lx 
-lcosh(0) 



(3) 



where 1 is the 2x2 unit matrix, f-i is the second Pauli 
matrix and 9 and x are position and energy dependent 
functions. We assume colinear magnetizations along the 
z-axis and s-wave singlet superconducting state. We 
choose a gauge where the superconducting order param- 
eter A is real and positive, and then the supercurrent is 
contained in the electromagnetic vector potential A and 
the chemical potential of the Cooper pairs is included 
in (j>. Inspection of the self-consistency relation for A 
reveals that \ — 0, it depending on the boundary con- 
ditions. This ansatz simplifies the calculations consider- 
ably. The advanced Green's function is related to the 
retarded through g A = — [p3g R P3] ^ ■ Because of normal- 
ization, the Keldysh Green's function can be expressed 
as g K = g R h — hg A where h is a diagonal distribution 
function matrix. 



We will now consider a stationary state. A kinetic 
equation can be derived from the Usadel equations if 
we include Keldysh components. The important quan- 
tities are the physical particle and energy currents (in- 
cluding particles and holes), which we will denote by 
jt and jl respectively, with the corresponding distribu- 
tion functions carrying the same indices, /it and /il [3]. 
The physical spin current is denoted jxs and the spin 
energy current jls, with distribution functions /its an d 
/ils- The spin-resolved distribution functions can be ex- 
pressed by the particle distribution function as /ixs(LS) = 
- (/T (E) ~ h (£))/2 - (+) (/ T (-E) - h i-E))/2. The cur- 
rent components jt etc. are spectral quantities, and the 
total charge current is given as an integral j c harge(r, t) = 
\e\No dEj^(r, t, E), and the spin current is obtained 
by a similar integral of jts- Energy current is given by 
joncr gy (r,t) = \e\N J^dE E} L (r,t, E), and the differ- 
ence in energy current carried by opposite spins by a 
similar integral of jls- 

The equilibrium solutions for the distribution functions 
are /il,o = tanh(/3i?/2) and /it,o = ^ls,o = ^ts,o = 0. 



We derive kinetic equations and find, 

V-j L = 0, (4a) 

V ■ j T = -2|A|a T T h T , (4b) 

V • j L s = - ^2|A|a T T + ^«lsls) /ils, (4c) 
1 
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(4d) 



The right-hand side terms represent rcnormalizcd scat- 
tering because of superconductivity: 



&tt =lm{sinh(0)} 



«lsls = (Re{cosh(fl)}r - (lm{sinh(0)}r 
<*tsts = (Re{cosh(6>)}) 2 + (Re {sinh(<2)}) 5 



(5a) 
(5b) 
(5c) 



The |A|axT terms on the right hand side in (4) are 
due to conversion of quasiparticle current into super- 
current, and the dLSLsAsf , aTSTs/ r sf terms are due 
to spin-flips. The spin- flip time in the normal state is 
r^ 1 = 87m s f7V S'(S'-|-l)|u3f| 2 /3, where n s f is the magnetic 
impurity density, Nq the density of states at the Fermi 
level, 5* the impurity spin quantum number and v s { is 
the Fourier transformed spin-flip impurity potential. We 
assume isotropic scattering. Our definition of r s f differs 
from the usual spin-flip lifetime by a renormalization fac- 
tor 4/3. This definition reproduces the diffusion equation 
with a spin-flip length l^' = \/Dt s [ in the normal state. 
Thus there is a difference between the spin-flip lifetime 
measured in e.g. electron spin resonance and spin-flip 
transport time. 

We introduce generalized energy-dependent diffusion 
coefficients 



Dt =D 



(Re{cosh((9)}) 2 - (Re {sinh(6<)}) s 



D T =D (Re{cosh(6»)}) 2 + (lm{sinh(0)}) : 



(6a) 
(6b) 



where D = TVp/3 is the diffusion constant. The currents 
can then be expressed as 

k = - D L Wh L +lm{j E }h T , (7a) 

j T = - D T Wh T + Im {j E } h u (7b) 

Jls = - -DtV/ils + Im {j E } /its, (7c) 

j'ts = - D L \7h TS + Im {j E } h LS . (7d) 

Here we have defined the spectral supercurrent as j_g = 
D{Vx ~~ 2eA) sinh 2 (0). The self-consistency relation is 



sgn(A )iVoA / d£sinh(0)/i L , (8) 



A(r) 



where the factor sgn(Ao) is determined from the bound- 
ary condition to give the correct sign and A is the inter- 
action parameter. The complex part of this equation is 
neglected as a consequence of charge conservation [8] . 
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The functions 9 and x are determined by the retarded 
components of the Usadcl equation. We obtain 



D 



\7 2 9 



-(Vx-2eA) 2 sinh(26») 



-2icosh((9)|A 



V-j B =0, (9) 
= - 2LEsinh(0) 
3 1 



4r s 



sf 



sinh(20), (10) 



where Eq. (9) implies that the spectral supercurrent is 
conserved. In addition we have the following symmetry 
conditions, 6*(-E) = -9(E), x*(~E) = x(E)- Equa- 
tions (4)-(10) determine all transport properties of S. 

In general, in a hybrid F/S system, the superconduc- 
tor cannot be described as in terms of BCS-formulas 
close to the F/S interface due to the proximity effect. 
Nevertheless, to gain insight into the physics implied by 
the above-mentioned formulas let us now consider the 
limit of a homogeneous BCS superconductor, and select 
X = 0. This is relevant for the proposed experiment 
below. For energies \E\ < |A| «tt = A/\/ A 2 - E 2 
and the spin-flip renormalization factors are cutsts = 
0, Q;lsls = -A 2 /(A 2 — E 2 ). The generalized diffusion 
constant Z?l = while D T = DA 2 /(A 2 -E 2 ). From Eq. 
(7a) this means that there is no energy current carried 
by quasiparticles with energy \E\ < |A|. Gap scatter- 
ing for quasiparticle energies below the superconducting 
gap corresponds to a transformation of the charge cur- 
rent (jx) into supercurrent. Such scattering is not possi- 
ble for the physical spin current (jxs)- Consequently, in 
the absence of spin-flip scattering the quasiparticle spin- 
current into the superconductor vanishes for \E\ < |A| 
since Z?l = in Eq. (7d) and Otsts = in the kinetic 
equation (4d). Note that this result relies on the fact 
that there are different effective diffusion coefficients for 
charge current (-Dt) and for spin current (-Dl)- We also 
observe that the term «lsls is negative below the gap, 
acting as a source of spin energy. 

Above the gap (\E\ > |A|) the factor 

aTT vanishes 

while a LSL s = £ 2 /(£ 2 -A 2 ), a TS TS = (E 2 + A 2 )/{E 2 - 
A 2 ). For the generalized diffusion coefficients we find 
that D L = D and D T = DE 2 /(E 2 - A 2 ). Now con- 
sider the kinetic equations in the BCS case. A charge 
current carried by quasiparticles with energy \E\ > \A\ 
can propagate into S. For quasiparticles at \E\ > |A| we 
see that there is no renormalization for the spin-energy 
diffusion length in Eq. (4c), whereas the spin diffusion 
length in Equation (4d) has an energy dependent renor- 
malization factor which diverges for energies \E\ = \A\ 
causing massive spin-flip scattering. 

We will now apply this formalism to study spin diffu- 
sion, and demonstrate the significance of the renormal- 
ization of the spin diffusion length. Experimental studies 
of spin accumulation and spin injection has recently been 
performed [9] in metallic spin valves. The spin accumu- 
lation in the physically different inelastic regime for a su- 



perconductor in this experimental system has also been 
calculated theoretically [7]. We will consider the simpli- 
fied geometry shown in Figure 1, where there is no charge 
transport in the superconductor, and calculate the spin 
accumulation signal in the elastic regime. The F1/N/F2 
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FIG. 1: Spin battery connected to a superconductor. The 
thick solid line indicates a tunnel barrier. 

systems act as a spin-battery which is connected via a 
tunnel junction to the superconductor. A voltage bias 
between Fi and F2 induces a spin accumulation that can 
flow into S. The superconducting wire is connected to an 
S reservoir in equilibrium BCS state by a good metal- 
lic contact at distance L^ s > from the N/S interface. On 
top of the S wire there is a ferromagnct connected by 
tunnel barrier which upon switching of the magnetiza- 
tion direction acts as a detector for the spin signal. Mea- 
surement of the relative voltage of this electrode between 
parallel and antiparallcl (with respect to the top F reser- 
voir) magnetization gives A/i = m' P ' — /^ Ap ) which de- 
scribes the difference between electrochemical potential 
of spin-up and spin-down quasiparticles located a dis- 
tance L from the N/S interface. This quantity can be 
calculated A/i = - dEP^h TS (L, E), where P( D) 
is the spin polarization of the tunnel barrier between S 
and the F detector. We assume a homogeneous order pa- 
rameter and BCS spectral properties in the S wire since 
there arc tunnel barriers between the N,F and S elements 
and perturbation from current and spin-flip is weak. 

We can express the difference between the spin-up 
and spin-down distribution functions in N close to S as 

A/ (N) = f W _ f (N) = p(FN) (/( £ _ e y /2) _ f ( E + 

eV/2)), where p( FN ) 



G^min) / (Cmaj ~t~ G m i n ) 



is the spin polarization between the F reservoirs and N, 
f(E ± eV/2) is the Fermi-Dirac distributions in the F 
reservoirs and G ma j( m in) is the conductance of major- 
ity (minority) spin electrons from ferromagnetic reser- 
voir to the middle of N. There is thus no charge current 
or supercurrent anywhere in S, however there may be a 
spin-current. Equation (7d) states that there is no spin- 
current for energies below the gap, thus for these ener- 
gies the N/S interface is effectively insulating. Since the 
S wire is connected to a reservoir in the other end for 
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\E\ < A the spin distribution function equals the equi- 
librium value /its = 0. We solve the TS kinetic equa- 
tion (4d) for energies \E\ > A. This equation reduces 
to a diffusion equation with renormalized spin-flip length 

4 S) (£) = W(£ 2 - A 2 )/(P 2 + A 2 ), where l sl = y/D^ 
is the normal state spin-flip length. The boundary con- 
dition at the S reservoir is that the distribution function 
attains the equilibrium value, and at the S /N interface we 
match at each energy the tunnel spin current to the spin 
current inside S, |e|./VojTS- We assume that £< s > /4 S) > 1 
which is a relevant physical situation. 

The position and energy dependent solution /its is sub- 
stituted into the expression for the measured difference 
in electrochemical potential for parallel and antiparallcl 
configuration, and we obtain 



/>oo 

A A t = 2p( D > / dPA/^c 

J A 



R 



(S) 



R 



(S) 



(11) 



where R m (E) = 1/(\T\ 2 N BCS (E)N ) is the resistance of 
the N/S tunnel barrier, \T\ is the tunneling matrix ele- 
ment, N BCS (E) is the BCS density of states, R { ^ } (E) = 
l^(E)p/A is the resistance of the S wire within a spin- 
flip length and p is the resistivity of the material in S 
when in the normal state (T > T c ). This result can be 
understood as follows. The spin accumulation close to 
the tunnel interface is exponentially attenuated by spin- 
flip scattering in S. The spin signal is also decreased by 
the tunnel resistance, and since spin current is strongly 
suppressed for energies \E\ < A only quasiparticles with 
energies higher than the gap contribute. The effective to- 
tal resistance is a series of the tunnel interface resistance 
with the resistance of S within one spin-flip length. 

We will now consider some simplified limits for 
the quantity A/i defined above. In the normal 
state where A — > we find that Afi/eV = 

2 p(D)p(FN )i? (S) cxp( __ L/Zsf)/(i? (S) + whcrc R (S) 

and assume their normal state (energy indepen- 

dent) values. At k B T <C A the signal measured 
by AyU vanishes when the bias is lower than the en- 
ergy gap eV < A since spin current is suppressed 
for quasiparticles below the gap. For higher bias, 
eV > A, and at zero temperature when the bulk resis- 
tance dominates, R^ 3> R^\ an approximate solution 
is A/x = 2P( D > p( FN >A c- L /' rf {e - Lr2 / 2 ^ /r - e"^ 2 ^ + 
^Lj2U(exi[r^/Lj2k{] - ctf[^L/2Q)}, where r = 
2A/eV. In this case the relation between the energy 
gap and the bias determines the magnitude of the spin 
signal, and the exponential decrease of the signal. 

The temperature dependence of A/i in the general 
case is given by decrease from a constant value above 
T c as the temperature approaches zero. An example 
of this behavior is shown in Figure 2. Here we have 
used the approximate temperature dependence A = 
1.76T c tanh(1.74 x /T c /T- 1). Our calculations show 
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FIG. 2: Temperature dependence of Ap/eV. We use R^' = 
in the normal state. For the dotted lines L/l s f = 6, 
and for the solid lines L/l et = 7. The bias eV is 0.1A(T = 
0), 3A(T = 0), 10A(T = 0) for the lower curve to the higher 
curve, respectively. 



that the spin signal decreases due to superconducting cor- 
relations. For a large energy gap the spin accumulation 
vanishes completely at low temperatures. These effects 
can be explained by suppressed subgap spin-current and 
massive spin-flip at energies close to the gap because of 
the superconducting correlations. 

In conclusion, we have presented a formalism to de- 
scribe elastic spin transport in superconductors with 
spin-flip scattering. We find different effective diffusion 
coefficients for charge- and spin-current. The spin- flip 
length is renormalized in the superconducting case, and 
at energies close to the gap there is massive spin-flip. As 
an illustration we compute the difference in electrochem- 
ical potential due to spin accumulation in an experiment 
sensitive to the renormalization of spin-flip length. 
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